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Abstract 

This short review is devoted to the problem of the equihbrium of stehar dynamical 
systems in the context of the Vlasov-Poisson model. In a first part we will review 
some classical problems posed by the application of the Vlasov-Poisson model to 
the astrophysical systems like globular clusters or galaxies. In a second part we will 
recall some recent numerical results which may give us some some quantitative hints 
about the equilibrium state associated to those systems. 
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1 Introduction 



Globular clusters and galaxies are concentration of stars whose physical char- 
acteristics are such that allow astrophysicists to model them as being in some 
equilibrium state for not too long time scales. The large number of stars and 
some properties of the gravitational interaction — which evidently play a key 
role for their dynamics — made the Vlasov-Poisson model a good candidate 
for their modeling. In this short review we present some results and problems 
of this field of research. 



2 The Vlasov-Poisson system in the context of stellar dynamics 



The Vlasov-Poisson model is a mean field approximation of the gravitational 
potential generated by a large assembly of point masses in the context of the 
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non dissipative kinetic theory. In order to apply this model to self-gravitating 
systems like globular clusters or galaxies, we have to pose three classical hy- 
potheses : First of all, it suffices to consider only Newtonian gravity for the 
mean field interaction between stars 0- Second, all stars are statistically equiv- 
alent, particularly they do have the same mass and they do not evolve. This is 
a drastic approximation but it could be relevant in a mean sense and during 
quite long evolution stages of these objects. Third, the considered systems 
must be non-dissipative on the time scales of interest. Nevertheless, it is well 
known that stellar dynamical systems could be dissipative mainly by two pro- 
cesses : On the one hand , such systems could contain gas which dissipates 
their total energy by dynamical friction. This feature must be considered for 
modeling the visible part of spiral galaxies. On the other hand it is well known 
that deviations from the mean-field forces due to the discreteness of the actual 
stellar distribution will alter the distribution of kinetic versus potential ener- 
gies as computed from the Vlasov-Poisson equations. In fact, this is a problem 
of time scales. The dynamical time Td of self-gravitating systems depends only 
on its mean density p via the relation ^ 1/ ^/Gp where G is the Newton 
gravitational constant (see (0) for details). It is typically the time taken by 
a test star to cross the system. Another duration is under interest : the time 
Tj. required for the mean velocity of the system to change by of order itself. 
This time has been estimated by Chandrasekhar for uniform systems to be 
proportional to the dynamical one via the relation Tr ^ N Td/ IniV (see (0) 
p. 189 for details). It is typically the time for which the dissipative process like 
encounters plays a key role in the system's dynamics. Galaxies are composed 
of at least = 10^ stars, therefore this time limitation is not relevant in this 
case. The situation is not so clear for globular clusters. The number of their 
components ranges from 10^ to 10^. If initial stages (during a few hundred 
of dynamical times ...) could be non dissipative and therefore potentially de- 
scribed by Vlasov equation, during late stages — for billions years clusters — 
encounters dissipate energy and stellar dynamicists introduce Fokker-Planck 
formalism. 



3 Equilibrium as steady state solutions of the Vlasov-Poisson sys- 
tem 

Self-gravitating coUisionless systems could be modelized by a phase space dis- 
tribution function / (r, p, t) and a mean field potential ip (r, t). Vectors r and 
p are the usual conjugated position and momentum which are elements of R"^ 

^ The general theory of relativity could also be considered in the context of the 
cosmological principle. The system takes then into account the scale factor a (t) of 
the Universe. This formalism is often used in the context of formation of the large 
structures of the Universe. 
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for d dimensional systems. The functions / and ip are coupled by the Vlasov- 
Poisson system, which is for d = 3: 



dl 
dt 



+{E, f}=Q 



where 



E = h rmlj 

2m 



P(r,t) 



Gm 
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Atp = AttG m j fir, p', t) dp' 
ip bounded at oo 



r — r' 



The quantity E denotes the mean field energy of a test star and plays a central 
role in this system. The function p (r, t) represent the mean mass density 
distribution of the system. Writing the Vlasov equation as above — i.e. using 
the Poisson brackets — makes trivial the well known result that every positive 
and normed function /<, of the mean field energy is a steady state solution 
of this system. A natural question, initially posed by S. Chandrasekhar in 
the middle of the last century, could then be : What are the properties of 
physical systems whose distribution function writes fo = fo (E) ? Although 
the answer of the inverse problem was well known by astrophysicists 0, the 
natural direct one was solved only fifty years later (jij), mainly by using in 
this context a difficult, but classical, mathematical theorem by Gidas, Ni and 
Nirenberg (jl|). The main ingredients of this result could be sketched as follow 
: If the distribution function depends only on the mean field energy, then the 
mean mass density depends on the position r only through the mean field 
potential 



In this case, Poisson equation writes Aipo = c po (i'o) where c is a positive con- 
stant. The physical context allows additional hypotheses : Newtonian gravity 
imposes that ipo (i") is a negative function bounded at infinity. The classical 
continuous limit let us consider that po (r) is a positive and continuous func- 
tion. The GNN theorem (jl|) then allows us to claim that tpo = V^odr]) and 
therefore the corresponding system has spherical symmetry in the spatial part 
of the phase space. Properties of the system in the velocity part of the phase 
space are more evident : The dispersion velocity tensor is clearly proportional 
to unity, therefore the system is said isotropic in velocity space. Spherical and 
isotropic steady state solutions of Vlasov- Poisson systems was intensively stud- 
ied in the context of stellar dynamics as classical models for globular clusters 

^ As a matter of fact, a spherical self-gravitating system must be associated to a 
radial gravitationnal potential which produces naturally a radial density profile by 
Poisson equation 



fo = f (E) ^ Po{r) = m J f + mipo j dp' = po (ipo) 
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and galaxies. Stability of such systems was proved for monotonic distributions 
functions in the linear case (see (j2|) for a whole presentation of this topic) and 
also in some general non linear cases (j^. 



Considering other potential isolating integrals of motion in the gravitational 
field, we can define other classes of steady state solutions of Vlasov-Poisson 
systems : This classical property of Poisson brackets is sometimes called Jeans 
Theorem in stellar dynamics (see (0) for instance). 



When the distribution function depend additionally on the mean angular mo- 
mentum modulus fo = fo{E,L'^), an extension of the GNN theorem shows 
that if the mass density is monotonic, then the system is always spherical. 
However, the tangential velocity dependence of = r'^vf makes the system 
anisotropic in the velocity space. The asymmetry of such systems is associated 
to an imbalance in the ratio of radial over tangential star orbits velocity distri- 
bution. Since the end of the 80's, stellar dynamicists have understood that such 
anisotropy could be at the origin of an interesting instability. As a matter of 
fact, if anisotropic spherical stellar systems are generally stable against radial 
perturbations , it could be proved that systems whose cannot be infinitesi- 
mally perturbed by radial disturbances are intrinsically unstable(4;). This is 
the fine mechanism of Radial Orbit Instability : For pure radial orbit system, 
each star orbit extension is exactly the radius of the whole system. Thus, 
such a system cannot receive infinitesimal radial perturbation which affects 
by definition only an infinitesimal part of the system. On the contrary, any 
non radial perturbation associated to a given spatial direction could stretch 
or compress infinitesimally the spatial extension of an associated star orbit. 
This feature generates a tidal friction which makes an instability to grow and 
forms a triaxial system from an initial sphere. Radial Orbit Instability triggers 
when a sufficient amount of radial orbits is present in the system — a general 
criterion is given in (j3) from distribution function susceptibility to receive ra- 
dial perturbations. As indicated by numerical analysis (llll ). this feature could 
be at the origin of triaxiality in some self-gravitating systems like elliptical 
galaxies. 



Less is known about steady state solutions characterized by distribution func- 
tions depending, in addition of E, on more complicated integrals. If some 
special models associated to fo = f {E, L^) are clearly triaxial (see (0) for de- 
tails) , there is, up to now, no extensions of the famous GNN theorem which 
allows to claim anything in a general way. 
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4 Equilibrium of a stellar system 



Taking into account the time problem limitation presented in section [2] and 
their generally quiet physical properties, one can modelize globular clusters 
and at least elliptical galaxies by steady state solutions of the Vlasov-Poisson 
system. The fundamental question that we have to answer is : What are the 
associated distribution functions ? 

This problem was attacked by stellar dynamicists using three approaches : 
Thermodynamics, comparison with observational data and numerical simula- 
tion. 



The thermodynamical approach, e.g. (jl3l ). is based on the classical assumption 



of statistical physics which associates the equilibrium state to the maximum 
Boltzmann entropy one. In the late sixties, (js]) reconsider this problem and 
failed thoroughly : The classical result of these works is the isothermal sphere 
which distribution function is /o = / (E) oc exp {(3E). This failure comes from 
the fact that in the 3 — D case, such entropy maximizer no exists. Nevertheless, 
interesting features come from such an analysis when the system is put in an 
unphysical box — see (fiol ) for a review of this topic. Another problem is the 
fact that the Boltzmann entropy 



S = - J f In f drdp 



which is extremalized in this approach, is a conserved Casimir functional in 
the Vlasov-Poisson context [1] ! 

Observational approach consists in the integration in the models of observed 
properties of globular clusters or galaxies. A non exhaustive list of such a 
models is presented in classical text books like (jsl) or more recently (j?!). Most 
famous ones are the King model for globular clusters which is an arbitrarily 
truncated isothermal sphere, and the Navarro-Frenk-White(i8|) radial profile 
for dark matter halos associated to the galaxies. From a theoretical point of 



^ A Casimir functional is on the form C[f] = f C{f)drdp where C is a smooth 
function. It is well known that such functionals are time conserved quantities if / 
is a solution of the Vlasov equation: 



dt J df dt J df \ m dx 5x dp 

J I'm ox ox op J 

where the last equahty follows by intergrating the first term over x and the second 
term over p, since / — > as |x|, |p| oo. 
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view, all these models are poorly justified, except perhaps the Henon isochrone 
model (0) which is unfortunately generally unrealistic ... 

The last way is numerical. Taking into account the recent ability given to stel- 
lar dynamicist by numerical tools to produce realistic formation experiences, 
important advances were made at the turning of the last century. These anal- 
ysis are presented in the next section. 



5 A numerical approach to equilibrium properties 



The idea is to produce an equilibrium state from the gravitational collapse of 
an arbitrary physical set of point masses. This is not a recent idea and it was 
pioneered by Spitzer since the late sixties. The effective realization of such a 
study in a general way was the result of a serie of works displayed over more 



than twenty years ((jlll) and reference therein, (1121 ) and reference therein). In 
those works , a general understanding of observed equilibrium properties of 
globular clusters and dark matter halos of galaxies follows from a detailed 
analysis of numerical gravitational collapses of a sufficiently large set of N 
point masses. Gravitational collapse from physical initial state produces gen- 
erally a sphere. Flatness is possible by Radial Orbit Instability but it requires 
strong and inhomogeneous collapse q|. The radial mass density profile of such 
spheres splits into two distinct classes : Homogeneous initial conditions col- 
lapses into a core halo-structure. It consists of a constant density region (the 
core) which extension can reach the half-mass radius of the system. This core 
is surrounded by a radial power law decreasing density region (the halo). The 
other class is composed by sufficiently inhomogeneous initial conditions whose 
collapse toward a monotonic radial power law decreasing density without no- 
table central core. The origin of the collapsed core of such systems could be 
explained by the effect of the Antonov Instability discovered in the thermo- 
dynamical context (see ([3) for a review). These two classes could be directly 
linked to the two classes of Newtonian self-gravitating systems which are glob- 
ular clusters and galaxies. The classical hierarchical galaxy formation scenario 
is evidently linked to the inhomogeneous class. It produces collapsed core — 
which is a favorable process to form the observed and always mysterious super 
massive black holes — and potential flattening in violent cases. The smaller, 
quiet, more isolated and then homogeneous case, could be naturally inter- 
preted as the generic globular cluster formation process. This could explain 
in the same operation, their generic spherical shape, their typical core-halo 
density proflle, and flnally their generic lack of intermediate mass black hole 



An important amount of solid rotation could also be invoked but this property is 
generally not observed in a sufficient way in stellar systems, like elliptical galaxies 
in particular. 
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in central regioiH]- The case of spiral galaxies formation and evolution is more 
specific and the role of dissipative process cannot be neglegted. Perhaps Vlasov 
equation is no more relevant for the description of such structures. 



6 conclusion 



Taking into account some physical constraints, the globular clusters and the 
galaxies could be suitably modelized by the Vlasov- Poisson model. In addi- 
tion to this modelization, accurate numerical simulations allows to obtain a 
global understanding of important stages of their evolutions and of their main 
differences. However, a lot of fundamental works are always to be done before 
the study of very particular properties of the self-gravitating systems. 
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